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A Variational Principle for Dissipative Fluid Dynamics 
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School of Fundamental Science & Technology, Keio University, 
Yokohama 223-8522, Japan 

In the variational principle leading to the Euler equation for a perfect fluid, we can 
use the method of undetermined multiplier for holonomic constraints representing mass 
conservation and adiabatic condition. For a dissipative fluid, the latter condition is replaced 
by the constraint specifying how to dissipate. Noting that this constraint is nonholonomic, 
we can derive the balance equation of momentum for viscous and viscoelastic fluids by using 
a single variational principle. We can also derive the associated Hamiltonian formulation by 
regarding the velocity field as the input in the framework of control theory. 

§1. Introduction 

In many areas of physics, variational principles help us choose convenient vari- 
ables in formulating problems. The least-action principle can lead to the equations 
of motion not only for microscopic dynamics but also those for nondissipative macro- 
scopic dynamics. For example, in a perfect fluid, the sum of the kinetic energy and 
internal energy is conserved, and we can write down the action easily. There are 
also the conservation laws for mass and entropy; these holonomic constraints can 
be incorporated into the action by means of the method of undetermined multipli- 
ers. Thus, the Euler equation can be derived from the stationary condition of the 
action ^ ^ ^ ^ ^ ^ ^ ^ ^ 

The Navier-Stokes equation is the equation of motion for the Newtonian viscous 
fluid, and represents the momentum balance with the dissipative force taken into 
account in the linear phenomenological law.'^''^ Without this force, the equation 
is reduced to the Euler equation. The dissipative force can be derived using another 
variational principle called the maximum dissipation principle, where Rayleigh's dis- 
sipation function is minimized .'^''^ Using this principle amounts to assuming the 
linear phenomenological law. A set of variables with respect to which Rayleigh's 
dissipation function is minimized is usually different from a set of variables with 
respect to which the action is stationary.l^''i3 

One way of unifying these two variational principles is called Onsager's varia- 
tional principle,'^^ which has been used frequently in formulating dissipative dynam- 
ics of complex fluids .'^■'^'^^"^"^ Using this principle, we can obtain the balance 
equation of momentum with the dissipative force taken into account in the linear 
phenomenological law. However, applying this principle to the Newtonian fluid in a 
straightforward way, we can obtain only a linearized version of Navier-Stokes equa- 
tion, i.e., the Stokes equation.—'-J* 

A viscous fluid is not locally adiabatic, and the entropy conservation law cannot 
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be assumed. We instead impose a constraint specifying how to dissipate. Noting 
that this constraint is nonholonomic, unhke in the previous studies, we can derive 
the balance equation of momentum in the framework of a single variational principle. 
As mentioned above, this framework for a viscous fluid originates from that for a 
perfect fluid, which we briefly review in ^ Our main results are shown in S}3l where 
we apply our procedure to viscous and viscoelastic fluids. In SjH for completeness, 
we derive the associated Hamiltonian formulation in the framework of the control 
theory. Our study is discussed and summarized in the last section. 



§2. Perfect fluid 



We briefly review Ref. _5) with slight modifications. Let us consider the dynamics 
of a perfect fluid in a fixed container from the initial time tinit to the final time tfin- 
Let V denote the region occupied by the fluid. The velocity field v satisfies no- 
penetration condition at the boundary dV, i.e., 

n-v = OondV , (2-1) 

where n is the unit normal vector directed outside on dV. Let us write p for mass 
per unit volume, and s for entropy per unit mass. The internal-energy density per 
unit mass, e, is a function of p and s because of the local equilibrium. The action is 
given by the integral of the Lagrangian density with respect to the time and space 
considered. Apart from the constraints mentioned below, the Lagrangian density is 
given by 

£{p,v,s) = p\lv^-e{p,s)] , (2-2) 



^2 

which is the difference between kinetic energy density and internal energy den- 
sity. The variables are not independent because of the constraints mentioned be- 
low. Writing T and p for the temperature and the pressure, respectively, we have 
de = —pdp~^ + Tds in the thermodynamics, which leads to 

where the subscripts s and p indicate variables fixed in the respective partial dif- 
ferentiations. We write r for the time in the Lagrangian coordinates although it is 
equivalent to the time t in the non-relativistic theory. The partial derivatives with 
respect to r (dr) and t {dt) imply the Lagrangian and Eulerian time-derivatives, 
respectively. 

In the Lagrangian description, we label a fluid particle with its initial position 
a = (01,02,03), and write X = (Xi, X2, X^) for its position at time r. The time 
derivative of X denotes the velocity fields v, i.e., 

drX = V . (2-4) 

The endpoints of the path line are fixed by 

6X{a, tinit) = 6X{a, tgn) = . (2-5) 
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The volume element in the Lagrangian coordinates can be given by the determinant 
of the Jacobian matrix, 

J(a,r) = ^^l^^^. (2-6) 

We assume that J has no singular points in the space and time considered and that 
J(a, tinit) is unity. The conservation law of mass is given by 

p(a, r) J(a, r) - pinM = , (2-7) 

where Pinit(ci) is the initial value of mass density. The adiabatic condition implies 
that the entropy of a fluid particle in the perfect fluid is conserved along a path line, 

s(a,T) - Sinit (a) = , (2-8) 

where Sinit(ct) is the initial value of entropy density per mass. Equations (j2-7p and 
()2-8p give holonomic constraints. Using undetermined multipliers, 7, K and A, we 
can define the action as 

/ ?r / d^a { J£(p, s, v) +7 • {drX - v) +K{pJ - pinit)+AJ{s - Sinit)} , (2-9) 



which is a functional of 7, K, A, p, s, v and X. Taking the variations of Eq. ()2-9p 
with respect to these trajectories, respectively, leads to Eqs. (|2-4p . ()2-7p . ()2-8p . 

K = -^v'^ + e + ^, (2-10) 
2 p 

A = pT, (2-11) 

7 = pJv (2-12) 



5 f /I 2 .A dJ 



and 

Roman indices run from 1 to 3, and repeated indices are summed up. The sur- 
face integral terms, appearing on the way of this calculation, vanish because of the 
boundary conditions, Eqs. (|2-ip and ()2-5p . Calculating by means of the cofactors 
yields 

^'^ , (2-14) 



d{dXi/daj) dX, 

while some algebra yields 



Thus, we can rewrite Eq. ()2-13p as 

d d { fl 
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Substituting Eqs. ()2-10p and ()2-12p into Eq. ()2-16p . we successfully obtain the Euler 
equation in the Lagrangian description, 

dvt dp 

We can replace dvi/dr by {dt + v ■ V)vi to obtain the Euler equation in the Eulerian 
description, which is given by Eq. (|2^ belowP® 

In the Eulerian description, a fluid particle is labeled by the Lagrangian co- 
ordinates A = (741,^2)^3)5 which depend on the spatial position x. Since the 
Lagrangian coordinates are conserved along the path line, we have 

_Ai = -V ■ VAi . (2-18) 

The endpoints of a path line are fixed by 

6A{x, W) = 6A{x, tfin) = , (2-19) 



instead of Eq. ()2-5p , as discussed in Ref . i5j) . The mass conservation law Eq. ^ 
can be rewritten as 

p{x, t) = J-\x, t)pinit{A{x, t)) , (2-20) 
where the inverse of J is given by 

J-.(.,,) = filli2li2l. (2.21) 

The action in the Lagrangian description, Eq. ()2-9p . is rewritten in the Eulerian 
description as 



/ dt d^x {C{p,s,v)+pi{dtAi + v-VA,) + K{p- p;^itJ~^) + A{s-Sirnt)] 
Jtinit Jv 



(2-22) 

which is a functional of /3, yl, s, v and A. The stationary condition of Eq. (|2-22p 
yields Eqs. ([ZHD, (f2T0D . (f2TTD . (f2T8D . (fMD . 

pv + l3iVAi = (2-23) 

and 

-A = - V . (ft.) - J a4- + 9i^r^-'*9(9l^/-^a4-- ^2-'') 

Here, pinit and Sinit depend on A, but p and s are independent variables. Using 
Eqs. dm), duo]), dMI]), 

9(9^i/axj) dAi ^ ' 

and 
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we can rewrite Eq. ()2-24p as 



d ^ „ , „ . ( dK ^ ds\ dxj 
-ft = -V . (ftt,) +{p—-pT, ^ ^ 



dt \ dxj dxj J dAi 

Dividing Eq. (|2-23|) by p, we calculate the Lie derivative of the resulting equation 
with respect to v to find the sum of the Lie derivative of v and 

{VA,)(^ + v.v] (2-28) 



dt J \ p 

to vanish, with the aid of Eq. ()2-18p . The Lie derivative of v was given in Appendix A 
of Ref . EI) • Substituting Eq. (|2-27|) into Eq. (j2-28p , we can derive the Euler equation, 

p\^v + \vv'^ -vx{V xv)\ = -Vp . (2-29) 



^dt 2 

§3. Dissipative fluids 

3.1. Simple dissipative system 

To show the essence of our procedure, we consider a damped harmonic oscillator 
in a heat bath. We define the position of the oscillator, g, so that it vanishes in the 
balance. We write m and k for the mass of the oscillator and the spring constant, 
respectively. The internal energy of the heat bath, E, is a function of its entropy, 5, 
and its temperature is given by T = dE/dS. We assume that the frictional force, 
/, is exerted on the oscillator. The action is given by the integral of the difference 
between the kinetic energy and the sum of the potential energy and the internal 
energy over the time interval [tinit , ifin] > 

+ E{S)^^ , (3-1) 

which is a functional of trajectories q{t) and S{t). The frictional force converts the 
mechanical energy into heat. Thus, the realized trajectories should satisfy 

considering that the heat bath is always at the equilibrium. Any trajectory making 
the action stationary under the constraint of 

f5q + T5S = (3-3) 

satisfies Eq. (|3-2p because the infinitesimal variations 5S and 5q can be those in 
the direction of the time. Change in the entropy generally should depend on how 
the oscillator is moved in the time interval from tiait to tgn, in terms of the ther- 
modynamics. Thus, for example, the constraint of (5g = at t^^ in the variational 
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calculation does not always mean that of 6S = at ign- Hence, then the constraint 
Eq. (|3-3p is generally nonholonomic because it cannot be rewritten in a form that the 
infinitesimal variation of a function, U, of q and S always vanishes, i.e., 5U = 
The variation of the action Eq. ()3-ip yields 



dt i ( -m^ -kq\6q- T6S } = . (34) 



Here, we fix both ends of the trajectory of i.e., 

<5g(timt) = Ht^n) = . (3-5) 
Substituting Eq. ()3-3p into Eq. ()34p gives the balance equation of momentum, 

= -kq + / . (3-6) 

The frictional force, /, is given by —^dq/dt with ^ being the friction coefficient, 
according to the linear phenomenological law. This simple model can be regarded 
as a simple Kelvin- Voigt model of the viscoelasticity. If the damping follows the 
Maxwell model, the second term on the left-hand side (Ihs) of Eq. (|3-'2p should be 
replaced by the frictional force multiplied by the deformation rate of the viscous part 
(dash-pot). This model is discussed in ^3.31 

Let us consider whether a straightforward way of applying the method of unde- 
termined multipliers to the constraint Eq. ()3-2p works well. Modifying the action so 
that it contains the product of the multiplier and the Ihs of Eq. (|3-2p , we encounter 



, {^d^q, fdT dT\ 







V 





+ [TSSfl:-^ = ' (3-7) 



where we assumed that T is a function of S, in calculating the stationary condition 
of the modified action. Because the last term on the Ihs of Eq. ()3-7p does not always 
vanish under Eq. ()3-5p . it is found that we cannot derive Eq. ()3-6p in this way. 



3.2. Viscous fluid 

For a viscous fluid, we write aij and Jq for the viscous stress tensor and the 
heat flux, respectively. Let us deflne the rate-of-strain tensor Cij as {djVi + diVj)/2, 
and the equation corresponding to Eq. ()3-2p turns out to be 

pT {dt + v- V)s - aijCij + V • Jq = . (3-8) 

We assume the no-slip boundary condition 

v = on dV , (3-9) 

instead of Eq. ()2Tp . and thus, we can assume no dissipation at the boundary. We 
assume that the whole fluid is enclosed by an adiabatic wall, i.e., that no heat flux 
cannot pass through the boundary. It enables us to integrate Eq. (|3-8p over V to 
obtain 
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Here, / = V • cj"^ is the dissipative force per unit volume, with the superscript 
indicating the transposition. In the Lagrangian description, as in Eq. (j3-3p . we thus 
impose the nonholonomic constraint, 



[d^a {J{pTSs + f-SX)} = , (3-11) 
Jv 



instead of the adiabatic condition Eq. (|2-8p for the perfect fluid. Thus, the action, 
Eq. ()2-9p . for the perfect fluid can be rewritten into the action for the viscous fluid 

as ^ 

/ 'dr / d^a {JCip, s,v)+j- (drX - v) + K{pJ - p,^,^)] , (3-12) 

which is a functional of 7, K, p, s, v and X. The stationary condition of this 
action under the constraint, Eq. (|3-lip . can be obtained, as Eq. (l3-6p is derived from 
Eqs. and ([231) • Thus, we obtain Eqs. (I^TUIl . IFZH\i and 



Substituting Eqs. (|2-lUp and (j2-12p into Eq. (|3-13p . we successfully obtain the balance 
equation of momentum for the viscous fluid in the Lagrangian description. 

As in ^ we can formulate the variational principle in the Eulerian description. 
Let us rewrite Eq. (j3-10p as 

d^x {pT{dts + v ■SIs) + f -v} = Q . (3-15) 
Jv 

From Eq. (|2-18p . we can derive 

dxj dAi 



dAi dt ' 



(3-16) 



which can also be derived from differentiating x(A(x,t),t) = x with respect to t. 
By means of Eq. (j3-16p . the nonholonomic constraint Eq. (j3-lip can be rewritten as 

^,=.{,„,-||(,T^ + /,)m.}=0. (3.17) 

The action, Eq. (j2-22p . can be rewritten as 

[ *^dt [ d^x {C{p, s, v) + /3i{dtAi + V ■ VAi) + K(p - Anit J"^)} , (3-18) 
J Unit Jv 

which is a functional of (3, K, p, s, v and A. Calculating the stationary condition of 
this action under the nonholonomic constraint, Eq. (|3-17p . we can obtain Eqs. (|2-10p . 
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([248]) . (ICTjl . (12^ and 



d ^ „ , ^ . ( OK ^ ds „ \ dxj 



9i V 5xj J dAi 

1 , 5p ^ ^ dXj 



= -V.(ft.) + ^-p-g^+g^-/,jg^. (3.19) 

Following the same procedure as shown just below Eq. ()2-27p . we find that Eqs. ()2-18p . 
(12^ and ([3091) yield 

d \ 1 

—V + - 2 - V X {V yiv) \ = - Vp + V • ct"^ , (3-20) 
(jij 2 J 

which is equivalent to Eq. ()3-14p . The viscous stress tensor aij is a symmetric tensor 
because the angular momentum is conserved, as shown in Appendix[Xl If the viscous 
stress tensor a is given in the linear phenomenological law, each of Eqs. p-14p and 
()3-20p represents the Navier-Stokes equation. 

3.3. Viscoelastic fluid 

We consider a viscoelastic fluid. We still write e^j for the rate-of-strain tensor of 
the fluid, while we write Eij for the strain tensor of the elastic part and aij for the 
viscous stress tensor of the viscous part. The internal energy, €e, is a function of p, 
s and Eij, and satisfles 

dee = —pdp~^ + Tds + p~^KijdEij , (3-21) 

whereby Kij is defined. We assume the rate-of-strain tensor of the viscous part to be 
given by 

eij = eij -{dt + v- V)Eij . (3-22) 

We can assume that Eq. (|3-8p holds if the second term on its Ihs is replaced by 
—aijCij. In the Eulerian description, the resultant equation is 

/ d^x {pT {dts + v-Vs) + {dkCTjk + ^kidjEki) Vj + a^jdtE.j} = . (3-23) 
Jv 

Using the same procedure in ^3.2\ we have the nonholonomic constraint. 



Ceip,v,s,Eij) = pl^v'^-eeip,s,Eij)\ . (3-25) 



The Lagrangian density is given by 

1 

2 

Replacing £. in Eq. ()3-18p by Ce, we obtain the action, the stationary condition of 
which with respect to 6Ai turns out to be 

d „ „ . ( dK ^ ds d&ik dEui\ dxj 



dt V dxj dxj dxk dxj J dAi 

1 dv^ dp dujk \ dxj 

2 dxj dxj dxk J dAi ' 



^ ^ ; I dv^ dp dajk\ dxn 
V • if3^v) + -P-— + -ii - -1^ — I , (3-26) 
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with the aid of Eq. ()3-24p . as Eq. ()2-27p is derived. With respect to 5Eij, we have 

l^ij ~ ^ij ■ (3-27) 

Using the same procedure in ^3.21 we have the balance equation of the momen- 
tum in the Eulerian description, which equals Eq. (|3-2Up with fjjj replaced by aij. 
Calculations in the Lagrangian description can also be performed as in ^3.21 

Within the linear phenomenological law, assuming aij to be linear with respect 
to Cij and Kij with respect to Eij, we can find 

{dt + v- V)a^j (3-28) 

to be linear with respect to etj — etj. Here, eij is linear with respect to 
This gives the viscoelasticity of the Maxwell typcl^^'E^ 

Equation (|3-22p represents how the strain tensor is convected. Under mathe- 
matical notation, the last term on the right-hand side (rhs) of Eq. ()3-2ip should be 
p^^K^dEj, for example. If Cij, eij and Eij are respectively replaced by e*, e*- and 
Ej in Eq. (|3-22p , the way of convection of the contravariant tensor involves the 
upper convected time derivative. Hence, then what we have considered here is the 
viscoelasticity of the upper convected Maxwell type within the linear phenomeno- 
logical law.l^ 

3.4. Two-component fluid 

We consider a two-component fluid. Let pa and pb be the mass densities of 
components a and 6, respectively. The total mass density is given hy p = pa + Pb, 
and the mass fraction of component a is ^ = Pa/p- We define the chemical potential 
as the Gibbs energy per unit mass, and write p for the difference given by subtracting 
the chemical potential of the component b from that of a. The internal energy per 
mass is a function of pa, Pb and s, for which we write €2{pa, Pb, s)- We have 

de2 = ^dp + pdilj + Tds , (3-29) 
p^ 

which can be rewritten in terms of dpa, dpb and ds with the aid of 

dp = dpa + dpb (3-30) 

and 

dip = dpa dpb ■ (3-31) 

P P 

The components a and h have their respective velocities Va and Vb and La- 
grangian coordinates Aa and Ab. The mean velocity of the fluid is given hy v = 
ipVa + (1 — tp)vb. The equations corresponding to Eqs. (|2-2ip and ()3-16p can be 
respectively written as 

j-i ^ d{Aci,Ac2,Ac3) 

" d{xi,X2,X3) 

dxcj dAci . , 



and 
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where c = a,b. The Lagrangian density is given by 



Pe2iPa,Pb,s) . 



(3-34) 



In general, the dissipation is associated with the fluid viscosity, the heat con- 
duction and the mutual diffusion. Here, we consider a polymer solution, where the 
components a and b are a polymer and a solvent, respectively!^ For simplicity, we 
neglect the polymer elasticity and fluid viscosity, although they can be considered 
as in the preceding subsections. We have 

pT {dt + V ■ V) s + ^ ■ {va - Vh) + V ■ Jg = , (3-35) 

where ^ denotes the frictional force. As we derive Eq. ()3-17p . we can obtain 



c=a,b \ J y j 



(3-36) 



where and ^f, are introduced to satisfy ^„ = — = ^. As in ^3.21 the action is a 
functional of (3^, Kc, Pc, s, Vc and Ac with c running from a to b, and given by 



tfin f 

dt / d^x 



^2 + ^ {^ciidtAci + Vc ■ VAci + KciPc 



Pc init 



c=a,b 



, (3-37) 



where pc init is the initial mass density of the component c. Under the nonholonomic 
constraint Eq. (|3-36p . we find the stationary condition of Eq. ()3-37p to be 



Pc — Pc initJc ^ — 
PcV + Pci'^Aci = 

d 



and 



Kr 



dKc 



Vc-v ] + e-\ h/ic l 

p \ p 



(3-38) 
(3-39) 

(3-40) 



(3-41) 



where fia and pb are introduced to satisfy pa = —pb = P- As in ^3.21 we can derive 
the balance equations of momentum for the components. 



Pc 



d 1 

—v + -Vv^-VcX (Vxv)) 



Pc^ PaPb^ . ^ 
P P 



(3-42) 



The sum of Eq. (|3-42p over c = a,b gives Eq. (ICTj) with V • a'^ deleted. The 
resultant equation describes the momentum balance of the solution as a whole. The 
calculations in this subsection are in the Eulerian description; those in the Lagrangian 
description can be performed as in 
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§4. Hamiltonian formulation 

Applying the control theorjfiSJ''2Sli to the dynamics of the perfect fluid, we derived 
its Hamiltonian formulation.'^ Similarly, with slight modifications as shown below, 
we can derive the Hamiltonian formulation for dissipative dynamics. In the control 
theory, a subsystem controlled by another subsystem is called a plant. Let q and u 
denote the plant state and the input to the plant, respectively. We assume that the 
plant dynamics follows 

^Q(t)-F(g(t),n(t)) = 0, (4-1) 

whereby a function F is defined. In addition to q and n, we assume that the La- 
grangian, L, also depends on a variable S*, which will denote the entropy in the later 
application to dissipative dynamics. In the control theory, we define an evaluation 
functional ^ 

rdtL{q,u,S), (4-2) 

and determine the optimal input u* so that the evaluation functional is minimized 
under the constraint given by Eq. (|4-ip . We also require Eq. (|3-c{p and T = —dL/dS. 

We take the constraints into account using the undetermined multiplier, p, which 
is called costate in the control theory. The functional to be minimized, which we 
call action, can be defined as 

!^L{q,u,S)+p(^^^-F{q,u)^^ . (4-3) 

Let u\q,p,S) denote the input minimizing Eq. ()4-3p on condition that q^p and S 
are fixed. This means that the time integral of H = —L(q, u, S) +pF(q, u) from tjnit 
to tfin is stationary with respect to u at u = u^{q,p, S). Let us define 

Hiq,p,S)^H{q,J{q,p,S),p,S) , (4-4) 

which is usually called the Hamiltonian. The preoptimized action is defined as 

''dt l-H{q,p,S)+p^a\ , (4-5) 



dt 



which is a functional of q, S and p. The optimal input u* is given by u'^{q,p,S), 
when q, p and S satisfy the stationary condition of Eq. ()4-5p under Eqs. ()3-3p and 
(|3-5p . We can usually assume d'^L/{dSdu) to vanish, and thus is a function of 
only q and p. In this case, we find 

dH dL ^ , 

Using Eqs. ()3-3p . ()3-5p . ()4-5p and ()4-6p . we obtain a set of modified Hamilton's 
equations, which are Eq. (j3-2p . 

dq dH 
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and 

In a straightforward way from Eqs. ()3-2p . ()4-7p and ()4-8p . we have dH/dt = 0. 
Further discussion is given in Appendix [XI 

Let us apply the formulation above to the viscous fluid discussed in ^3.21 In the 
Eulerian description, the state, costate, entropy density and input are respectively 
given by A, /3, s and v. Equation (j4-ip is given by (j2-18p in this fluid, where the 
function corresponding to H is given by 

HiA^v, p, s) = -p I ^v^ - e{p, s) I + ■ VA, . (4-9) 

Here, p is regarded as a function of A given by Eq. (j2-2Up without the undetermined 
multipher K introduced. Let v'^{A,(3) denote the input corresponding to u'^, and 
we find the Hamiltonian to be given by the sum of the kinetic and inertial-energy 
densities, 

n{A,l3,s)^p!^^v^^ + e{p,s)^ , (4-10) 

where satisfies Eq. ()2-23p . The modified Hamilton's equations can be found to be 
Eqs. ()2-18p and ()3-19p . In deriving these equations, we note that, for example, 

6(3i d^i dxj d{djPi) ^ ■ ' 

should replace d/dp in Eq. ()4-7p because the fields are involved here unlike in 
Eqs. (j4-7p and (|4-8p . In the Lagrangian description, v'^ satisfies Eq. (|2-12p . and 
the modified Hamilton's equations are found to be Eqs. ()2-4p and ()2-16p with v 
replaced hy vK 



§5. Discussion 



It is natural that the nondissipative equation can be derived by means of the 
least-action principle. We can formulate dissipative dynamics by adding the dissipa- 
tive force to the equation of motion to this nondissipative equation. In this formu- 
lation, the variational principle is inherent to the nondissipative dynamics. We can 
obtain the dissipative force by means of the maximum dissipation principle, where 
the linear phenomenological law is simply derived from the stationary condition of 
the quadratic form. These two variational principles have been utilized especially in 
formulating dynamics of complex fluids, as mentioned in 311 

In the irreversible thermodynamics for a one-component viscous fiuid,'^''^ we 
start with the mass conservation law, the balance equations of momentum and en- 
ergy. Then, combining these equations with the local equilibrium yields Eq. (j3-8p 
from which the rate of the total entropy production is written in terms of fiuxes and 
thermodynamic forces.^ The linearity between them is assumed to give the linear 
phenomenological law, determining the dissipative force near the equilibrium. 
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Suppose that we do not know the balance equation of momentum but know 
Eq. (j3-8p . representing how the dissipation occurs. Still, we can calculate the rate 
of the total entropy production from Eq. (|3-8|) . and then we can determine the 
dissipative force within the linear phenomenological law. On the other hand, we can 
derive the balance equation of momentum from Eq. (j3-8p by retracing a part of the 
calculations mentioned in the preceding paragraph. What we do in this work is to 
formulate this derivation in a single variational principle. 

This variational principle is easy to understand in terms of physics. The least- 
action principle for the perfect fluid is associated with the entropy conservation law. 
Replacing this constraint on the entropy with a constraint describing how the dissi- 
pation changes the entropy, we can obtain our variational principle. The key point 
inherent to the dissipative dynamics is to appreciate that this is the nonholonomic 
constraint to which the method of undetermined multipliers cannot be applied. 

Applications of this variational principle to a viscous fluid, a viscoelastic fluid and 
a two-component fluid are shown in ^ They are straightforward generalizations of 
the simple case discussed in ^3.11 We believe that this principle can be easily applied 
to the dynamics of other complex fluids, such as liquid crystals and electromagnetic 
fluids.El.ED 
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Appendix A 

Symmetry 

We show that the Hamiltonian version of Noether's theorenP^liiEHli g^jn works 
for the model discussed in ^ Let us consider a canonical transformation, map- 
ping (qojPoiSo) to {qa,Pa, Sa), which makes the modified Hamilton's equations, 
Eqs. ()3-2p . ()4-7p and ()4-8p . invariant. When the transformation with respect to 
a is infinitesimal, the generator G{qo,pa, Sq, t) can be defined. Then, because of the 
canonical transformation, we have 



/ dt {pogo - Ho{qo,po, Sq)} = dt < Paia - Ha{qa,Pa, Sa) + ^ {qoPa + OlG) \ , 

(A-l) 

where {qo,Po, Sq) and {qa,Pa, Sa) respectively satisfy the modified Hamilton's equa- 
tions. With the aid of Eqs. ([321) and ([331), Eq. (TATI) yields 



5q = qo 



qo 



a 



dG 

dpo 



(A-2) 
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A - f dG fdG\ 

and 

do 

Haiqa,Pa, Sa) - Ho{qo,po, So) = a— , (A-4) 

where we simply wrote G for G{qo,po, So,t). From Eqs. p-3p . (|4-6p - (j4-8p . ()A-2p and 
()A-3p . we have 

Haiqa,Pa, Sa) = Haiqo,Po, Sq) + -^^Q + ^"^^^ + ^^'^ ^^'^^ 

fdpdG dqdG dSdG\ 
= H^(^o,Po,So)-a^-— + -— + -—j . (A.6) 

Substituting Eq. (jA^ into Eq. yields 

HaiqoiPo, So) - Hoiqo,Po, So) = a— . (A-7) 

Thus, the condition dG/dt = is equivalent to the invariance of H, i.e., = Hq, 
under the canonical transformation. If time-invariant, H is the generator of the 
infinitesimal canonical transformation with respect to the time. 

As discussed in ^ the Hamiltonian density of a fluid is given by T-L{A^ (3, s) of 
Eq. (j4-10p . where A, (3 and s are respectively the state, costate and entropy density. 
As Eq. (]A-7p . we can derive 

jd^x {Ha{Ao,Po^so)-Ho{Ao,f3o,so)} = jjPx . (A-8) 

Thus, if Ti is invariant under the transformation generated by Q, the integral of Q 
over V is conserved. This means that a vector field J can be defined so that 

^+V-J = (A-9) 

is satisfied and that the surface integration of J-n over dV vanishes. Its converse also 
holds. The spatial transformation is generated by the momentum, pv^ = —f3iVAi. 
The total momentum of the fluid is conserved when there are no external forces 
acting on the boundary dV. It is thus natural that we can rewrite Eq. (|3-20p in the 
form of Eq. ()A-9p . i.e., 

l(pvl) + d,{pvlvl-Tij) = 0, (A-IO) 

where Tij = p5ij — dij is a stress tensor. The spatial rotation is generated by the 
angular momentum, I = px x = —x x /3jVAj. The total angular momentum of 
the fluid is conserved when there is no external torque acting on the boundary dV. 
The time evolution of I turns out to be 

^ + ^ {^i^! - eijkXjTk?! - eijkTkj = , (A-11) 

with the aid of Eq. (jA-lOp . where the Levi-Civita symbol eijk is anti-symmetric on 
each pair of indic es. T he last term on the Ihs above should vanish from the discussion 
just around Eq. ()A-9p . which leads to aij = aji. 
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